Phase-noise induced limitations on cooling and coherent evolution in opto-mechanical systems 
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We present a detailed theoretical discussion of the effects of ubiquitous laser noise on cooling and the co- 
herent dynamics in opto-mechanical systems. Phase fluctuations of the driving laser induce modulations of the 
linearized opto-mechanical coupling as well as a fluctuating force on the mirror due to variations of the mean 
cavity intensity. We first evaluate the influence of both effects on cavity cooling and find that for a small laser 
linewidth the dominant heating mechanism arises from intensity fluctuations. The resulting limit on the final 
occupation number scales linearly with the cavity intensity both under weak and strong coupling conditions. 
For the strong coupling regime, we also determine the effect of phase noise on the coherent transfer of single 
excitations between the cavity and the mechanical resonator and obtain a similar conclusion. Our results show 
that conditions for optical ground state cooling and coherent operations are experimentally feasible and thus 
laser phase noise does pose a challenge but not a stringent limitation for opto-mechanical systems. 
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Over the past years tremendous experimental progress with 
opto-mechanical devices OJ, 0] and analogous systems in 
the microwave regime JH |jl has been made. Many groups 
have by now achieved significant cooling of mechanical mo- 
tion Jl H, 0, H 0, Gil 111 and in several of these setups 
resolved-sideband conditions have been demonstrated, which 
is a necessary prerequisite for quantum ground state cool- 
ing |[l2i IT3I ll4L Hal . Beyond laser cooling, which works 
well in the regime of weak opto-mechanical coupling, the re- 
cently demonstrated strong coupling regime 1 16| might allow 
for the observation of coherent dynamics between the cavity 
field and the vibrating mirror. For example, the transfer of 
single photonic excitations to the phonon mode has been sug- 
gested [TTtI [Titl to prepare and study quantum superpositions 
of macroscopic objects. 

Despite a steady experimental progress cooling to the 
ground state and further the combination of strong cou- 
pling conditions with low occupation numbers have not been 
achieved so far. The main limitation in current systems arises 
from the re-thermalization rate of the mirror Y m - k^T /hQ,„, 
(where Q,„ is the mechanical quality factor and T the tem- 
perature of the support), which competes with optical cool- 
ing. However, with lower base temperatures and increas- 
ing mechanical quality factors, mechanical heating can be 
strongly reduced and the impact of other imperfections on 
opto-mechanical systems must be considered. In particular, it 
has been argued recently lfl9ll that ubiquitous laser phase fluc- 
tuations impose severe limitations on opto-mechanical cool- 
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ing schemes and can impede ground state cooling under cur- 
rent experimental parameters. While effects of laser noise 
have indeed been observed |0], the laser linewidth require- 
ments Ti * 10~ 3 Hz for ground state cooling estimated in lfl9ll 
are clearly inconsistent with current experimental achieve- 
ments where residual mean occupancies as low as «o - 30 
lEol I2H l22ll have been achieved with Tl ~ 1 kHz. This dis- 
crepancy can be understood from a suppression of noise at the 
mechanical frequency a) m |0], but a rigorous derivation has not 
been given so far. Furthermore, potential impairing effects of 
a finite laser linewidth on strongly coupled opto-mechanical 
systems have not been addressed yet. 

Given the great interest in opto-mechanical experiments the 
impact of phase noise in such systems deserves a thorough 
theoretical investigation, which we provide in this work. We 
here generalize the standard descriptions of opto-mechanical 
systems to include laser phase noise with arbitrary frequency 
noise spectrum S^(Q) and evaluate its influence on cooling 
and coherent oscillations. The main results of this work are as 
follows. In accordance with previous predictions |0 [T^] we 
obtain a lower limit for the final mirror occupation number 
which is proportional to the intensity and the noise spectrum 
S x{u) m ) at the mechanical resonance frequency a> m . Surpris- 
ingly, this result applies for both the weak and strong cou- 
pling regime and we find that for a given T m and go (opto- 
mechanical coupling per single photon ), the condition 

S^(0) m ) < gl/r m , (1) 

must be satisfied in order to achieve ground state cooling. In 
the strong coupling regime we derive a similar condition for 
the observation of coherent oscillations, 

Sj,(aj m ) < gl/K, (2) 
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where k is the cavity field decay rate. Importantly, we also 
show that under relevant conditions decoherence from low 
frequency phase noise in opto-mechanical systems is negli- 
gible. Thus a small laser linewidth Ti/k <k 1 together with 
suppression of phase noise at high frequencies is sufficient to 
enable ground state cooling and the observation of coherent 
oscillations. The conditions (Q]) and (f2]i are experimentally 
challenging, but well within reach with state of the art laser 
stabilization. 

In the next section we develop the model for opto- 
mechanical coupling including laser noise, closely following 
the lines of IU3l 12311 . The results on weak and strong cou- 
pling are derived in Sec. [n] and Sec. [Ill] respectively. Details 
of calculations are moved to three appendices. 

I. MODEL 

We consider a typical opto-mechanical setup of an optical 
cavity mode coupled to a micro- or nanomechanical oscillator. 
We will refer in the following to the prototype example of an 
opto-mechanical system, a Fabry-Perot cavity of length L with 
one fixed, heavy mirror (which serves as the input coupler) 
and a light vibrating micro-mirror. However, our results apply 
as well to microtoroidal cavities or to dielectric membranes 
in a Fabry-Perot cavity. We restrict our analysis to a single 
vibrational mode of the mirror of mass m and frequency a> m 
and a single cavity mode of frequency co c that is driven close 
to resonance by a laser of frequency a>i. According to lf24ll 
the field couples to the motion of the mirror via the radiation 
pressure interaction and a total Hamiltonian can be written (in 
a frame rotating at a>f) 

H = -A a* a+co m tfb+goa?a(b+b f )-i (£* (i)a- fi(/)a f ) . (3) 

Here a, a* are the bosonic operators for the cavity mode 
such that [a, a'] = 1 and b, £>' for the mechanical mode 
[b, b^] = 1, while Ao = 0)i - u> c . The optomechanical cou- 
pling is go = (u> c /L) y/h/ (2mu> m ) and £(f) is the noisy laser 
field with an absolute value |£(f)l ~ So- While our analysis 
can be generalized to arbitrary noise processes we focus in 
the following solely on phase noise, 

&{t) = 8 eW . (4) 

Here 0(f) is a Gaussian noise process with zero mean which is 
characterized by the correlation function 

{^)^')) c/ = J ^(Q )e -^-'>, (5) 

where {0} ct denotes the average over different noise realiza- 
tions. The frequency noise spectrum S ^(£2) is specific for each 
experimental setup. For concreteness we will consider below 
a simplified noise model with a noise spectrum and correlation 
function, 

W = J^a- N^')) d = r L7ce -y^\ (6 ) 



respectively. Here F^ can be identified with the laser linewidth 
which is typically known. The additional parameter y" 1 char- 
acterizes a finite correlation time of the underlying noise pro- 
cess and leads to a suppression of noise at high frequencies. 
In the limit y c — > oo we obtain the <5-correlated noise model 
used in 01]. 

In addition to the classical driving field, the cavity field is 
coupled to the electromagnetic modes of the environment. We 
describe the resulting dissipative dynamics by a master equa- 
tion 

p = -i[H,p\+£ f (p), (7) 

where -C/(p) = K(2apa^ — a"' ap- pa* a) and k is the cavity field 
decay rate. We made use of the fact that for optical experi- 
ments even at room temperature the thermal photon number 
is effectively zero lE^ll . In principle Eq. (0 should also con- 
tain the term describing coupling of the mechanical system to 
its thermal environment. The influence of thermal noise on 
opto-mechanical cooling schemes has been studied in previ- 
ous works lfl2l [l3l [Til [l5tl . To focus on the effects of laser 
noise only we here consider the regime of efficient laser cool- 
ing where thermal noise is already suppressed and the dynam- 
ics is well described by the master equation (0. 

Linearized opto-mechanical coupling 

The radiation pressure from the cavity field leads to a mean 
displacement of the mirror and we are only interested in fluc- 
tuations around this shifted equilibrium position. We there- 
fore perform a unitary transformation b — » /3 + b where 
fi = -gon p h/(o m is the mean displacement amplitude and 
yiph = {{a' i a)q}ci is the mean cavity photon number averaged 
over both quantum and classical fluctuations. In addition, we 
perform another unitary displacement operation for the intra- 
cavity field a — > a{t) + a where 

a(t)= f dse- { - iA+K)(, - s> 8(s), (8) 

is the classical part of the stochastically driven cavity field and 
A = Ao - 2go/3 is the effective detuning. Note that by choosing 
this time-dependent displacement amplitude a(t) instead of a 
constant one |[l9ll we can unambiguously interpret a as the 
quantum field of a non-driven cavity. 

From these transformations we obtain the master equation 

p = -i[H opt (t) + g Q N(t)(b + b\p] + £ f (p), (9) 

where N(t) = \a(t)\ 2 - {laWPlrf measures intensity fluctua- 
tions of the mean intracavity field and the opto-mechanical is 
now given by 

H ovt (t) =u m b*b - Aa'^a + g (a*{t)a + or(?)a t ) (b + />"'") . (10) 

Note that in Eq. ( fTOb we have already neglected the nonlinear 
term ~ a* a which is equivalent to an expansion in 1 / and 
well justified by current experimental setups. 
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From the equation of motion (O we conclude that laser 
noise contributes two effects: First, it causes a modulation 
of the opto-mechanical coupling strength ~ a(t). Second, it 
induces a stochastic force on the mirror ~ M(t), which arises 
from a conversion of phase to amplitude fluctuations inside 
the optical resonator. This heating mechanism has no direct 
analog in laser cooling in free space but a similar effect has 
been observed with trapped atoms inside a cavity ll26il27ll . In 
the following we evaluate the influence of both contributions 
on cavity cooling and the coherent evolution in this system. 



the stochastic quantity e"* (/) , which itself is a non-linear func- 
tion of the noise process (f>{t). For the remainder of the paper 
we therefore assume that the total change in phase tp{t) ac- 
cumulated on a timescale of k is small. This assumption is 
equivalent to the experimentally relevant limit F^ <K k and 
allows a rigorous expansion of noise correlation functions in 
powers Tl/k. The details of these calculations are shifted to 
App . lAl and IB] and here we summarize only the main results. 

Amplitude fluctuations 



II. WEAK COUPLING REGIME 

We first consider the weak coupling limit go\a(t)\ <K k. 
In this regime the cavity mode acts as a dissipative chan- 
nel for the resonator mode which for example can be ex- 
ploited for cooling. We are interested in the effects of phase 
noise on the mean mirror occupation number in steady state 
«o = {(b'b) q } c i (t — > oo). From the linearized Hamiltonian 
([Tol l we can derive a set of coupled differential equations for 
the quantum averages, e.g. (b^b) q , which can be solved in 
the perturbative limit go\a(t)[_«: k. Averaging the result over 
the classical noise process 12811 we obtain an effective cooling 
equation 



« = -W(n - «o) , 



(11) 



with n = {(b*b) q } c i. Here the cooling rate W = S(u) m ) - 
S(-U) m ) and the average steady state occupation number 
«o = S(-a> m )/W depend on the total noise spectrum S(a>) 
S n(oj) + Sa(w) which is the sum of the modulated intracavity 
amplitude correlation spectrum 



o 



dr {a*(T)a(0)} cl e iAT e-" T e 



(12) 



and the spectrum of intensity fluctuations, 

/-*oo 

S N (cS) = 2gfae dr {N(t)N(0)} c1 e m . (13) 
Jo 

For the case of a non-fluctuating driving field S(t) = So, 
a(T) = ao = &o/(iA + k) and S n(oS) = 0. Then, by evaluating 
Eq. (fT2l for the optimized detuning A = A op = - Jk 2 + w 2 , 
we recover the standard results, 



Wo 



2goVo| 2 



|A 



1 /I VI , 

no = 1 

2 \ (Jj m 



(14) 



which are well known from the theory of opto-mechanical 
laser cooling ITLZl [l3"ll . In particular in the sideband re- 
solved regime a> m » k the minimal occupation number is 
no = k 2 /4ojI. 

Our goal in the following is to study the influence of a noisy 
driving field on the final temperature. From the definition of 
a(t) given in Eq. ([8j and Sa(m), S n(cj) given in Eqs. ( fT2l 
and ([T"3l we see that the final occupation number « depends 
on integrals over two and four point correlation functions of 



Let us first look at phase noise induced modifications of 
Sa(w)- To provide some physical insight we can assume in a 
zeroth order approximation that the cavity field simply follows 
the driving field adiabatically, 



a(t)*S e m /( K -iA). 



(15) 



By inserting this expression into the definition of Sa(cl>) we 
obtain 



dT{e'^ T) e-' m 

o 



(16) 



and a simple result for this integral can be found in the white 
noise limit y c — > oo where {e'* (T) e _ '^ (0) } c ; - e~ ri ' T L From this 
estimate we expect correction to 5a(±W;h) and therefore to the 
final occupation number «o which are of order 0{Ti/k). 

In App.lAlwe present a more rigorous calculation of Sa{<jJ) 
for Fl/k <sz 1 and here briefly discuss results for the resolved 
sideband regime k <k <d,„. For transitions on the red sideband, 
which are associated with a resonant exchange of a vibrational 
quanta and a cavity photon, we obtain 



Wq 



~ 1 



CdQ ' 
J 2n k 



+ o 2 



> 1 - — . (17) 



Here we obtain a dependence on the low frequency part of the 
noise spectrum and the lower bound which has been derived 
using the model defined in Eq. (|6]l agrees with the estimates 
from above. For the blue sideband transitions, which corre- 
spond to a simultaneous excitations of a photon and a phonon 
we obtain instead 



Wo 



K 

4a& 



18^2 



(2co m ) 



(18) 



i.e., a modification which depends on phase noise at high fre- 
quencies. The contribution at 2o) m can be understood form the 
fact that for A = —Oi m it takes this amount of energy to excite 
a photon and a phonon. Although noise at to m is non-resonant 
with this excitation process it is maximally enhanced by the 
cavity response, and therefore leads to an equivalent contribu- 
tion. However, in both cases the noise is either non-resonately 
driving the system or is suppressed by the cavity response and 
the resulting corrections to the final occupation number only 
scale as ~ /c 2 /w„ xTl/k, i.e. they are reduced by the sideband 
parameter k/oj,,, . In summary we conclude that phase modula- 
tions of the cavity field lead to a reduction of the cooling rate 
by Fl/k, the ratio of laser to cavity linewidth, and therefore 
pose no serious experimental limitation. 
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Intensity fluctuations 

We now look at intensity fluctuations. As the zeroth order 
approximation of a(t) would lead to a constant intensity |or(f)l 2 
we expect the first correction for a slow noise process to scale 
as ~ 4>(t)lK and we therefore approximate 



a(t) « — 1 1 - 1- 



k- iA 



k- ;A 



(19) 



Inserting this expression into the definition of Sn{<jJ) in 
Eq. (f]~3T > we obtain for A = -a> m , 



Wo 



f 

Jo 



a | 2 Re {4>(T)4>{<d)}cie^" T dT. (20) 



Already form this simple estimate we find that in accordance 
with earlier predictions Jgl [Ht] intensity fluctuations do im- 
pose a limit on the final occupation number which increases 
with the mean cavity photon number l^ol 2 - However, Eq. (1201 
also predicts a crucial dependence of this limit on the partic- 
ularities of the noise process, i. e., the correlation function 
{4>{t)4>{Q)} c i. This means that the final occupation number de- 
pends in a non-universal way on the phase noise characteris- 
tics and cannot be inferred from a white noise model described 
by the linewidth only. 

The basic prediction of Eq. ( f20b is confirmed by a more 
rigorous derivation outlined in App.|B] where we evaluate the 
correlation function \N(t)N{Qi)} cI to first order in the param- 
eter Fl/k. Under this assumption and A = A op the resulting 
limit on the final occupation number can be written as 



no > 



W~ 



I for 



,S^(w ffl ) /|A, 



2k 



(21) 



This is the main result on laser noise induced limitations of 
opto-mechanical cooling in the weak coupling regime. It 
agrees with a simplified analysis given in Ref. f9] for the 
resolved sideband regime and generalizes the result for the 
white noise limit derived in Ref. Il9tl for arbitrary noise pro- 
cesses. It is instructive to consider the toy model for laser 
noise as given in Eq. (|6]l from which we obtain in the side- 
band resolved regime 



with the known bounds from the theory of opto-mechanical 
cooling with an ideal laser lfl2l[l3ll we get 



2kT„ 



"o 



2 + 4ZJ + 



(23) 



Here T m = k B T/tiQ,„ is the mechanical heating rate for a res- 
onator with quality factor Q,„ coupled to a thermal phonon 
reservoir of temperature T. The first term on the right hand 
side is the residual thermal occupation, the middle term is the 
contribution from heating due to Stokes scattering, and the 
last term stems form laser noise. For an optimal choice of 
the intracavity photon number |aol 2 = (2K) 2 r m / gfa ^(uim) we 
obtain 



"o 



\F m Sx{(jJ m ) k 



2 

Aiol 



(24) 



From this result we identify Eq. ([]]) as the relevant condi- 
tion to achieve ground state cooling in weakly coupled opto- 
mechanical systems, provided the optimal value of |ao| 2 is not 
prohibitive. Note that this condition depends on go, the radia- 
tion pressure coupling per single photon, cf. Eq. (O. 

The details of the laser noise characteristics will depend on 
the concrete experimental setup, but the spectrum of intensity 
fluctuations S w(w) can in each case be directly measured, e.g., 
from correlations of the transmitted intensity 7 O ut(0- In the 
strongly driven regime the cavity out-field is dominated by the 
classical part b out (t) ~ s/Kajt) and using standard results on 
photon counting statistics 12911 we obtain the simple relation 



Wo 



Si(a>) 



- 1 



(25) 



^drCii^e'^ is the spectrum of the 
normalized photon current correlation function C/(t) = 
{/out(T)/ O u,(0)W{/out(0)} 2 - 1 and S m is the shot noise con- 
tribution thereof. 



III. STRONG COUPLING REGIME 



n > \a \ 



2 L 



Tc 



« 7c + W 2 



(22) 



For typical experimental parameters \ao\ ~ 10 , a>,„ ~ 10 
MHz, k x 1 MHz and V L * 1 kHz the assumption of a 
white noise model, y c — > oo, would lead to the prediction 
«o Z 10 5 , which is in sharp contrast with observed experi- 
mental data. However, for a realistic noise model with a finite 
cutoff frequency y c «: lo,„ intensity fluctuations are strongly 
suppressed and the resulting limits no > 1-100 are consistent 
with recent experiments. 

It was pointed out already in ijgll that the scaling of the 
bound in (fJTJ with the intracavity photon number |»ol 2 implies 
an optimal driving power balancing the cooling effects with 
heating due to laser noise. If we combine the present result 



We now consider an opto-mechanical system operated in 
the strong coupling regime where the linear photon-phonon 
interaction strength G = goao exceeds the cavity linewidth 
k. This regime has been recently studied experimentally 
lfl6ll and has been discussed theoretically in the context of 
opto-mechanical cooling ll23i l30lll . However, more impor- 
tantly strong coupling conditions enable a coherent exchange 
of photonic and mechanical excitations. Thus, the opto- 
mechanical system can serve as a quantum interface between 
photons and phonons with potential applications for state 
preparation and quantum measurements of macroscopic me- 
chanical motion. It is therefore worthwhile to study the role 
of phase noise in particular under strong coupling conditions 
where due to the required large values of |ao| more pro- 
nounced effects are expected. 
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Strong coupling 

Let us first briefly review the main features of the opto- 
mechanical system in the strong coupling regime, ignoring 
for the moment the presence of phase noise or other imper- 
fections. When the interaction between the cavity and the res- 
onator mode exceeds the cavity decay rate the system dynam- 
ics is conveniently described in terms of the new collective 
operators A ± which diagonalize H opt , 



H opt = w+A+A + + w_AlA_ . 



(26) 



Assuming resonance conditions A e = -a> m and |G| < co m /2 
which is required for stability fTill . the eigenfrequencies are 
w± = o) m (l + 2\G\/a> m ) 1 ^ 2 and the normal modes are approxi- 
mately given by 



» m G* le" , 
: a H — — b 



y[2 \/2 2w ffl \ yfe 



— b'±—a ] '\, (27) 



V2 



where we have defined e' 2e := ao/|ffol- For n °t too large val- 
ues of \G\ the eigenmodes are essentially equal superpositions 
of the original cavity and resonator mode and they are split 
in frequency by w + - w_ 2\G\. A normal mode splitting 
exceeding the cavity decay rate k is a first signature of the 
strong cou plin g regime and has recently been observed in ex- 
periments Il6ll . In this limit the Liouville operator is approxi- 
mately given by l23ll 



-Cfip) + l) J] { 2A (P A l - A \ A & -p a I a z) 



+ ^"0 Yj [ 2A \P A t - A t A \P - P A t A \) ■ 



(28) 



As expected we see that both modes decay with half of the 
cavity decay rate. A lower limit on the achievable occupation 
numbers «o = \G\ 2 /4u> 2 n <K 1 arises from small admixtures 
of a ' and b T in Eq. d27l ) due to energy non-conserving terms 
in H opt . Neglecting this small correction we immediately see 
from Eq. ( f28b that without any additional heating mechanisms 
the total number of excitations in the system n to[ = (A|A + ) + 
(ALA_) decays as 



«,nt = -KH„ 



(29) 



Therefore, in the strong coupling limit the opto-mechanical 
cooling rate W is independent of the driving strength and sat- 
urates at the maximum value set by the cavity field decay rate 

K. 



Strong coupling in the presence of phase noise 

In the presence of noise the picture above is modified on 
one hand by the presence of fluctuating forces and on the other 
hand by a modulated cavity -resonator coupling G — » G(t) : = 
goa(t). To account for this time-dependent coupling we intro- 
duce time-dependent mode operators A ± (f) which diagonalize 



the Hamiltonian H opt (t) given in Eq. ( [Tol l at each point in time, 
i.e. 



[H p t (/),A ± (/)] = -co ± (t)A ± (t). 



(30) 



Here a> ± (t) = u> m (l ± 2\G(t)\/(o m )^ 2 are the instantaneous 
eigenfrequencies and the decomposition of A ± (f) contains 
now also time dependent phases e' 2e(,) := a(t)/\a(t)\. To sim- 
plify the following discussion and to focus on the effect of 
noise we assume that the condition k «: \G(t)\ <k lo,„ is strictly 
fulfilled, which allows us to neglect in Eq. d27b contributions 
of order 0{G{t)/a> m ). Under this assumption 



p = -i[H(t),p]+£ f (f)(p), 



(31) 



where 



(32) 



and 



■Cf(t)ip) * \ Yu ( 2A ?( f )P^(f) - A|(*)A f (*)p - P Al(i)Af(ij) . 



(33) 

We now investigate how those modifications affect cooling 
and coherent dynamics in the strong coupling regime. 



Cooling 

We first look at phase noise induced limitations for ground 
state cooling in the strong coupling regime. In the previous 
discussion on the weak coupling regime we have seen that 
the effect of the modulation of a(t) is less crucial for cooling 
than that of the fluctuations of the cavity intensity. We there- 
fore neglect for the moment the explicit time dependence of 
A ± (f) and study the effect of JV(f) only. A justification for this 
approximation follows form the analysis present below. We 
obtain the coupled equations 



(N ± ) = -K(N ± )-ig N{t)(e- m A 
<A ± > = -(icj ± + K /2){A ± )-ig N(t)e 



e W A ± ), 



(34) 
(35) 



where A^ ± - A ± A ± . After averaging over the noise we end up 
with 



{(N ± )} c i = -k{(N ± )} c , + S n (oj ± + iK/2), 



(36) 



where S a<(w) is defined in Eq. dT3l . In the limit Fl/k « 1 we 
can use the results derived in App.|B]to evaluate this quantity. 
By assuming that for frequencies D. ~ <o m the noise spectrum 
5 is flat on a scale k and to lowest order in k/G we obtain 



S N (a> ± +iK/2) ^ g \a a \ 



4ar m S j,(co ± ) S j,(u) m ) 



(tuT ± - Ojl) 2 2(<jl) ± - (L) m ) 2 



(37) 



Interestingly, heating arises from two contributions. As ex- 
pected, the first term represents intensity fluctuations at the 
eigenfrequency a> ± . Since in the strong coupling regime these 
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frequencies are well separated from the cavity resonance, the 
noise at these frequencies is suppressed by u) 2 J{<x> 2 ± - co 2 n ) 2 - 
The second contribution in Eq. (l37l i represents noise which is 
resonantly enhanced by the cavity. Although it is not resonant 
with co+ it still couples to the damped motion of A ± and we 
see that both terms lead to similar contributions for heating. 
For oj ± oj m ± G the lower limit for the total steady state 
occupation number is then given by 

«tot > ml ■ (38) 

K 

Surprisingly for S^(co m ) ~ S ^(cj±) this result is quite similar 
to the weak coupling regime which we attribute to a cancella- 
tion of two effects. On one hand for increasing org the cooling 
rate saturates in the strong coupling regime at W — K. One 
would therefore naively expect a scaling no ~ |aol 4 - However, 
since the eigenfrequencies a) ± are well detuned from the cav- 
ity resonance the effect of phase fluctuations is suppressed by 
1/G 2 ~ l/gg|aol 2 ( see Eq. ( |37] i) which reduces in the strong 
coupling limit the scaling from |aol 4 to |ao| 2 . 

Due to the saturation of the cooling rate the strong cou- 
pling regime does not offer a particular advantages for cool- 
ing. We therefore do not go further into details and study in- 
stead the effects of phase noise for coherent dynamics where 
strong coupling conditions are essential. 

Coherent oscillations 

As already mentioned above, the interesting aspect about 
the strong coupling regime of opto-mechanical systems is the 
ability to realize a coherent interface between mechanical and 
optical modes. For example, if the resonator is initially pre- 
pared in the ground state (possibly in combination with other 
cooling methods) and the cavity in a Fock state |l) c this sin- 
gle excitation is swapped onto the resonator mode at time 
t s - n/(2\G\) and recovered at a later time t r = n/\G\. More 
general we can describe this process in terms of an arbitrary 
coherent state \% c ) c which in a frame rotating with ai r evolves 
under ideal conditions as 

|0>,|&> e -» I - ie- i2 %) r \0) c -» |0) r | - % c ) c . (39) 

Since this evolution is independent of the initial coherent state 
amplitude it can be generalized to arbitrary quantum states. 
Therefore, apart from a known phase the dynamics generated 
by H op implements a faithful mapping between the states of 
the cavity and the resonator mode. If this operation is fast 
compared to photon loss, i.e. k <k \G\, it can be employed as 
a coherent way for an optical preparation of motional states. 
Alternatively the reverse process would enable an optical de- 
tection of the resonator state. 

To characterize coherent oscillations in the presence of pho- 
ton loss and phase noise it is sufficient to study the evolution 
of an initial coherent state |^o) = \&(Q))r\€c(Q))c under the 
evolution of the effective Hamiltonian H e jf(f) = H(t) - iKC^c. 
Ignoring small corrections of order G / <±> m <k 1 which are not 
essential in the following discussion the state will then evolve 



into \if/(f)) = \^r(t)) r \^ c (t)) c where the corresponding coherent 
state amplitudes can be written as 

(&(t)\ _ (cbbit) c ba (t)\ (&(0)\ (c b (t)\ ,, m 

~ \Cab(t) C an (t)j \U0)J + \c a (t)j ■ {W ' 

The ideal evolution described in Eq. ( f39b suggest to use 
the state overlaps |<-jf c (0)c -afi |f c (0>| 2 and |<-£.(0)|£.(f)>| 2 to 
characterize a state transfer or a full oscillation respectively. 
We here choose the latter option and for an distribution P(r, 9) 
of initial coherent state amplitudes £.(0) = re"* we define the 
fidelity 

f{i) = -\ rdr dd>{e- m+re '* l \iP(r,G). (41) 
x Jo Jo 

To be more concrete we average over initial states with r — 1 
and for £ r (0) = 0we finally end up with 

rit) = e -{\c m mn 2 h xe -t\cM\\,_ (42) 

While there is some arbitrariness in this choice of a fidelity, 
the definition (l42l) is sensitive to different aspects of the noise 
and should therefore be a good characterization for coherent 
processes involving a low number of excitations. For a noise- 
less system c a {t) — and we obtain 

T (t) = e -li+-s(|c|,) C -»/ ¥ _ (43) 

The fidelity for a full oscillation is approximately given by 
To(t r ) ^e- 7 * 2 * 2 ^ 2 . 

To study the additional degrading of T(f) in the presence 
of noise we write a(t) =* e~ m (A + (t) - A_(/))/ V2. The coef- 
ficients c aa {t) and c a {t) can then be calculated from the evolu- 
tion of the mode operators, 

.(A + \ (oj + (t)-iK/2 -6»(f) \(A + \ Kf ^( e ~ m \ 

(44) 

From this expression we identify three potential effects of 
phase noise. First, intensity fluctuations ~ N(t) introduce a 
random displacement c a (t) which is related to the additional 
heating discussed above. Second, the fluctuating cavity am- 
plitude \a(t)\ leads to a fluctuating normal mode splitting, 
u) ± (t) ^ oj ± ± Saj(t) where 6cj(t) = goN(t)/\ao\. Third, phase 
fluctuations of the intra-cavity field causes non-adiabatic tran- 
sitions between the modes A ± proportional to 6{t). 



High frequency noise 

Based on our discussions so far we expect that intensity 
fluctuations are the dominant decoherence mechanism. In a 
first approximation we therefore neglect the time dependence 
of co{t) and 6{t) and study the effect of N{t) only. By integrat- 
ing Eq. ( 1441 we then obtain 

r(t) =f X e- D( '\ (45) 



7 



where 



D(t) = {\c a (t)\\, 



-k(t+t')/2 



~-2gl f dr f dT'e m '" (T - T,) e- K{: 
Jo Jo 

x sin(|G|T) sm(\G\T'){N(T)N(T')} c 



(46) 



For a full oscillation period t = tt/(\G\) and |G| » k this ex- 
pression reduces to 



D(t = n/(\G\)) * 8|or 



(47) 



As a consequence the error e = 1 - f for a full oscillation 
between cavity and resonator mode is 



golaol 



+ 8|o | 



(48) 



For fixed cavity parameters there is an optimal field ampli- 
tude |q-()| for which, apart from a numerical prefactor, the error 
scales as 



fs 



(49) 



We see that apart from other imperfections, achieving 
S ${u) m ) <K g^/K is a necessary requirement for the observa- 
tion of coherent dynamics in opto-mechanical systems. Since 
k must also exceed the mechanical heating rate F,„ this result 
implies that the observation of coherent oscillations puts more 
stringent bounds on the acceptable level of phase noise than 
just cooling. 



Low frequency noise 

While due to the scaling ~ |ffo| 2 intensity fluctuations im- 
pose a sever limitation on coherent oscillations, this effect 
depends on phase noise at relatively high frequencies a> ~ 
u> m . Since the relevant system dynamics occurs on a slower 
timescale |G| this noise can in principle be filtered out in a 
carefully designed experimental setup. An important question 
therefore is whether or not other decoherence mechanisms ex- 
ist which depend on low frequency regime of the phase noise 
spectrum. 

To address this question we now assume that the phase 
noise spectrum S^(Q) has relevant contributions only at Q <K 
to m . This assumption allows us to omit the term ~ N(t) in 
Eq. ( PPfl i and to study the effects of 5a>(t) and 9(t) only. For 
a simplified discussion we will also neglect cross-correlations 
between these two stochastic quantities. In App. ICl we show 
that the resulting fidelity is then given by 



T(t) ~ exp(-{|l +cos(|G|f)e 



-<d/2 a -W(t)/2 e -R(t)/2 



/2 \ 2 )c)- (50) 



Compared to the case of a noiseless laser we obtain two addi- 
tional contributions to decoherence. Here 



W(t) = 2 



fds P 

Jo Jo 



ds' {6a)(s)6co(s')} c 



(51) 



describes dephasing due to a modulated normal mode splitting 
and 



R(t) = 2Re 



f ds f ds'e' 
Jo Jo 



G(i-s') 



{e(s)8(s')} cl , (52) 



describes non-adiabatic transitions between the normal 
modes. 

As anticipated above Eq. ( BTT l and Eq. d52b differ qualita- 
tively form the heating discussed above in the sense that the 
quantities W(t) and R(t) depend on noise far below the res- 
onator frequency. To see this more explicitly we use 5u>(t) = 
goN(t)/\ao\ and obtain 



W(t) 



goVol 2 



C dQ. 4sin 2 (Q 
J|n|< w „ 2*" ^ 2 



(«r/2) 



SAO). (53) 



An upper bound for this integral can be obtained form the long 
time limit, 



|G| 2 

W(t) < l -^S^(Q) x t. 

cot 



(54) 



We see a dependence on zero frequency noise S^O) ~ Tl, but 
since this noise is far detuned from the cavity resonance this 
dephasing process is suppressed by |G| 2 /w 2 . For the parame- 
ter regime of interest the resulting error for a coherent oscil- 
lation e x |G|r £ /w 2 «: 1 is therefore always smaller than the 
error due to photon loss. 

A similar conclusion can be obtained for decoherence 
caused by non-adiabatic transitions. A simple estimate for the 
upper bound on R(t) can be obtain by assuming that the in- 
tracavity phase 20(f) follows the laser phase tp(t) adiabatically 
(see Eq. (fT5]l). Under this approximation we obtain 



R(t)< is^(|G|)xr. 



(55) 



Again we see that for Ti <k k this decoherence process is 
negligible compared to cavity loss and a more accurate cal- 
culation predicts a further reduction due to a suppression of 
low frequency noise by the cavity response. We conclude 
that while low frequency noise does affect the dynamics of 
the opto-mechanical system it is for the parameter regime of 
interest negligible compared to photon loss. 



IV. SUMMARY & CONCLUSIONS 

In summary we have analyzed the effect of laser phase 
noise on cooling and coherent dynamics in a generic opto- 
mechanical system comprised of a driven optical cavity with a 
vibrating end-mirror. Our approach significantly extends and 
generalizes previous treatments to this problem and provides 
a rigorous way to study different aspect of the laser noise in 
such systems. For opto-mechanical cooling in the weak cou- 
pling regime our predictions our predictions for a final occu- 
pation number «o ~ O(10) are consistent with experimental 
data and we derive a condition for the noise spectrum which 
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is required in future experiments to achieve ground state cool- 
ing. We have also show that the discrepancy between exper- 
iments and previous theoretical calculations lfl9ll is based on 
the assumption of a white noise model, which does not lead to 
physically meaningful predictions. For opto-mechanical sys- 
tems in the strong coupling regime we have shown that the 
effects of phase noise still scale linearly with the intensity but 
the observation of coherent oscillations places more stringent 
bounds on tolerable level of phase noise. Nevertheless, we 
conclude that ground state cooling and coherent state transfer 
experiments can be achieved with state of the art laser stabi- 
lization techniques. 
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APPENDIX A: EVALUATION OF THE AMPLITUDE 
FLUCTUATION SPECTRUM 

We evaluate the amplitude fluctuation spectrum 5^(w) de- 
fined in Eq. dT2b in the limit Tl/k «; 1. From the definition 
of a{t) given in Eq. ^ we find that under stationary condi- 
tions the correlation function of the classical cavity field can 
be written in the form 

/~»oo c2y 

{a*(T)a(0)} c i = I dye' 2 ** I dx e~ iAx C 2 (T - x) , (Al) 

Jo J-2y 

where C 2 (t) = {£*(f)£(0)} c / is the two point correlation func- 
tion of the driving field. For Gaussian phase noise we obtain 



C 2 (t) = |£o|V? 



(0 2 W)rf 



ft 12 

where (p(t) = J_ t , 2 ds <f>(s). Using the noise model defined in 
Eq. © we see that the expectation value in the exponent is 
bound by {(f> 2 (t)} d < r L y c t 2 for t < y; 1 and {(f> 2 (t)} d < Y L t for 
long times. Therefore, in the limit Ti «c k we can expand the 
exponential in Eq. (IA2t to first order, 

/ i w/2 w/2 \ 

C 2 {t) - |£ | 2 (l - - J ds J ds' {i>(s)i>(s')} cl j . (A3) 

Equivalently, we can rewrite this expression in terms of the 
noise spectrum, 



C 2 (t) - |fio 



'(-/ 



dQ. 1 - cos(Or) 



(A4) 



We insert this result back into the definition of {a*(T)a(Q)} c i 
and Sa(<±>) and evaluate the remaining integrals. Since the 
resulting general expressions are lengthy we here only present 
the results for the sideband resolved regime k <k co m , where 
A op = -cj„, and Wo = 2g 2 ) \ao\ 2 /*"■ For the red sideband we 
obtain 



Wo 



r da ' 

J 2tt k 



Q 2 



(A5) 



Using the noise model defined in Eq. (O we see that in the 
white noise limit y c » k corrections are of the order of 
S ^(0) /(2k) * Tl/k while in the opposite limit they scale as 
TiJcIk 2 . For the blue sideband transitions we obtain the re- 
sult presented in Eq. ( fT8l . 



APPENDIX B: EVALUATION OF THE INTENSITY 
FLUCTUATION SPECTRUM 

We evaluate the intensity fluctuation correlation function 
Cm(t) - {N(t)N{0)} c i in the limit F^ «: k. Ignoring small 
corrections of order O(go/i<) the average cavity photon num- 
ber is n ph = {(a*a) q } c i {|ar(/)| 2 } d and N(t) = \a(t)\ 2 - 
{|ar(/)| 2 } c z- Under stationary conditions Cjv(t) can be written 
as 

C u {r)= dy x dy 2 e ' AV 
Jo 

-»2vi /*2y 2 
1-2 



X 



f yi dxi f dx 2 e~ iAx ' e~ iAx * C 4 (r-(y l -y 2 ),x u x 2 ). 

J-2y, J -2 



(Bl) 



Here we have introduced the four point field correlation func- 
tion 

C A (T,t u t 2 ) ={£* (T + h/2)&(T - h/2)S* te/2)£(-/ 2 /2)} e z 
- {£* (h 12) 6 (-h /2)}{8* fa/2) 8 (-t 2 /2)} d . 

(B2) 

For phase noise with Gaussian statistics the four point corre- 
lation function is given by 



C 4 (T,t u t 2 ) = \& \ 4 (e-i 



(<S 2 (7\fi,f 2 )) c , _ „-|{0 2 ft)+«S 2 (f 2 )) ( 



'). (B3) 



(A2) where <D(t, t , , t 2 ) = JT*f ds if>{s) + £' 2 ' 2 ds j>(s) and 0(r) 



+t 2 /2 



J_^ 2 ds<p(s). Following the same argumentation as in the 
evaluation of Sa(co) in App. [A] we can in the limit IV, «: k 
expand the exponentials in Eq. (IB3b to first order, 



C A {T,t u t 2 ) * -|£ | 4 ds \ ds'{j>(s)i>(s')} cl , (B4) 

Jr-t, J-t 2 

or in terms of the noise spectrum S^fi), 

sin(Qf j /2) sin(Q?2/2) _ iClT 



C 4 (T,t u t 2 )^^\& 



y J Tn 



D 2 



(B5) 

After inserting this expression back into Eq. (IB lb and evalu- 
ating the remaining integrals we obtain 



|4 fdn 4A%(Q) g -«* 
C N {r) - lorol j 2^ (A 4 + 2 AV-Q2) + («2 + 



Q 2 ) 2 )' 



(B6) 

For the evaluation of S N(a> m ) defined in Eq. ( TT3b the integral 
over t results in a term ~ S(Cl - oj m ) and we end up with 



Siv(w m ) = gl\ao\ 4 



4A 2 SJco m ) 



(A 4 + 2A 2 (k 2 - w 2 ) + (k 2 + w 2 ) 2 ) ' 



(B7) 
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For the optimal detuning A = - y/ajft, + k 2 this expression sim- 
plifies to S N ((o m ) = gl\ao\ A Sj,{co m )lK 2 . 

In the strong coupling regime heating rates depend on the 
fluctuation spectrum evaluated at the imaginary frequency, 
S ff(o) ± + ik/2). In this case we obtain 



2i„. |4 



S N (u> ± + ik/2) =g \a 



dQ. 



4k 



2n k 2 +4(Q- w ± ) 2 

4A 2 

X (A 4 + 2A 2 (k 2 - Q 2 ) + (k 2 + Q 2 ) 2 ) ' 

(B8) 

For A = -a> m and k <zz G this integral is dominated by two 
contributions from the resonances at Q. — a> m and Q. — a> ± . In 
the limit k — > we obtain the result presented in Eq. < T37T >. 



APPENDIX C: LOW FREQUENCY NOISE 

In the absence of intensity fluctuations, N{t) = 0, the for- 
mal solution of Eq. ( l44b is given by 

(i:l)-- iW (^- itM(i¥i )(l:(o))' < ci > 

where T is the time ordering operator, 



Mn 



CJ+ — ik/2 
cj- — ik/2 



(C2) 



and 



To obtain {c aa (t)} c i we take the classical average of Eq. ( IC1I ) 
and evaluate the average of the exponential of M(f) using a 
second order cumulant expansion. Then 



[A + (t)}d 
{A-(t)}c 



,-mt -m(o M+(0) 



U_(0) 



(C4) 



where 



M(/) = 



1 (W(t) + R(t) + il(t) -X*(t) 



X(t) 



W(t) + R(t) - il(t)j ■ (C5) 



Here W(t) and R(t) are defined in Eq. (P5TJ, and Eq. d52l ), and 
7(f) is the imaginary part of integral in Eq. d52l i which leads to 
a small shift of the oscillation frequency. Finally, 



X(t) = 2 f ds f ds'(e- w 
Jo Jo 



{Sco(s)6(s')} c , 



-i\C\ 



s {8a>(s')e(s)} cl ), 



(C6) 



is an additional cross term which arises from correlations be- 
tween 6co{t) and 9{t). However, using the a rough approxi- 
mation for a(t) » ao(l _ $(t)/(.K + i^m)) and 6{f) » 0(f) we 
obtain {5w(f)0(f'))d ~ \G\I Lo m {(f>{t)(f){t')} c i and also this term 
does not contribute significantly to decoherence. For a sim- 
plified discussion, both the cross term and the frequency shift 
are omitted in Eq. d50l ). 
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